Fast permutation preconditioning for fractional diffusion equations
نویسندگان
چکیده
In this paper, an implicit finite difference scheme with the shifted Grünwald formula, which is unconditionally stable, is used to discretize the fractional diffusion equations with constant diffusion coefficients. The coefficient matrix possesses the Toeplitz structure and the fast Toeplitz matrix-vector product can be utilized to reduce the computational complexity from [Formula: see text] to [Formula: see text], where N is the number of grid points. Two preconditioned iterative methods, named bi-conjugate gradient method for Toeplitz matrix and bi-conjugate residual method for Toeplitz matrix, are proposed to solve the relevant discretized systems. Finally, numerical experiments are reported to show the effectiveness of our preconditioners.
منابع مشابه
Solving large systems arising from fractional models by preconditioned methods
This study develops and analyzes preconditioned Krylov subspace methods to solve linear systems arising from discretization of the time-independent space-fractional models. First, we apply shifted Grunwald formulas to obtain a stable finite difference approximation to fractional advection-diffusion equations. Then, we employee two preconditioned iterative methods, namely, the preconditioned gen...
متن کاملPreconditioned Generalized Minimal Residual Method for Solving Fractional Advection-Diffusion Equation
Introduction Fractional differential equations (FDEs) have attracted much attention and have been widely used in the fields of finance, physics, image processing, and biology, etc. It is not always possible to find an analytical solution for such equations. The approximate solution or numerical scheme may be a good approach, particularly, the schemes in numerical linear algebra for solving ...
متن کاملA numerical method for solving a class of distributed order time-fractional diffusion partial differential equations according to Caputo-Prabhakar fractional derivative
In this paper, a time-fractional diffusion equation of distributed order including the Caputo-Prabhakar fractional derivative is studied. We use a numerical method based on the linear B-spline interpolation and finite difference method to study the solutions of these types of fractional equations. Finally, some numerical examples are presented for the performance and accuracy of the proposed nu...
متن کاملFinite integration method with RBFs for solving time-fractional convection-diffusion equation with variable coefficients
In this paper, a modification of finite integration method (FIM) is combined with the radial basis function (RBF) method to solve a time-fractional convection-diffusion equation with variable coefficients. The FIM transforms partial differential equations into integral equations and this creates some constants of integration. Unlike the usual FIM, the proposed method computes constants of integ...
متن کاملPreconditioning Techniques for Diagonal-times-Toeplitz Matrices in Fractional Diffusion Equations
The fractional diffusion equation is discretized by an implicit finite difference scheme with the shifted Grünwald formula, which is unconditionally stable. The coefficient matrix of the discretized linear system is equal to the sum of a scaled identity matrix and two diagonal-times-Toeplitz matrices. Standard circulant preconditioners may not work for such Toeplitz-like linear systems. The mai...
متن کامل